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OBJECTIVES 


1 .  Mraloata  the  currant  aataorological  measurement  technique  to  determine  if 
it  ia  sufficient  to  adequately  describe  the  surface  meteorological  processes 
used  to  infer  evaporation  duct  height. 

2.  Compare  relative  performance,  sensitivities  to  meteorological  inputs,  and 
ease  of  computation  for  the  evaporation  duct  height  models  developed  by  jeske, 
Mother am,  and  Davidson. 

3.  Compere  and  evaluate  the  EM  wave  propagation  path loss  models  developed  by 
Rotheram  and  the  Naval  Ocean  System  Center  (NOSC). 

4.  Determine  a  maximum  range  of  detection  error  for  the  modified  NOSC  propa¬ 
gation  model  as  employed  by  the  Integrated  Refractive  Effects  Prediction 
System  (XREPS). 


RESULTS 

1.  The  evaporation  duct  height  models  of  Jeske,  Rotheram,  and  Davidson  func¬ 
tion  equally  well  with  reasonable  duct  heights  calculated  when  considering  the 
statistical  average  of  meteorological  inputs. 

2.  The  meteorological  observations  taken  by  operational  personnel  are  not  of 
sufficient  quality  to  insure  proper  duct  height  calculations  when  individual 
duct  height  predictions  are  desired. 

3.  When  considering  statistical  means  in  observational  data,  the  pathloss 
models  of  Rotheram  and  NOSC  performed  reasonably  well. 

4.  The  assumption  of  a  homogeneous  atmosphere  for  a  fleet  operational  prop¬ 
agation  loss  model,  is  reasonable. 

5.  The  theory  behind  calculation  of  duct  height  and  pathloss  values  is 
statistically  sound. 


RECOMMENDATIONS 

1 .  Quality  of  individual  meteorological  observations  needed  as  input  to 
current  state-of-the-art  propagation  models  is  insufficient.  This  can  be 
improved  by  (a)  improving  the  observational  skills  of  the  individual  through 
further  training,  and  (b)  developing  a  sensor  system  that  can  be  deployed  away 
from  the  observation  platform. 

2.  Develop  an  entirely  new  duct  height  model  using  statistical  techniques. 

3.  Design  controlled  experiments  to  investigate  the  role  of  each  assumption 
upon  the  actual  pathloss  over  the  propagation  path. 

4.  Train  fleet  operators  of  the  IREPS  model  in  the  limitations  and  recog¬ 
nition  of  situations  in  which  the  model's  performance  is  marginal  or  poor. 

5.  undertake  further  evaluation  of  the  IREPS  model  by  comparing  a  predicted 
range  with  an  actual  observed  maximum  detection  range. 
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INTRODUCTION 


BACKGROUND 

The  ataosphere  is  in  a  constant  state  of  turbulence  because  of  differen¬ 
tial  solar  heating.  This  turbulence  produces  a  continual  variation  in  hori¬ 
zontal  and  vertical  temperature  and  aoisture  gradients. 

As  electrcsiagnetic  (EM)  waves  propagate  through  this  medium  of  varying 
temperature  and  moisture,  they  refract,  or  "bend."  A  measure  of  the  degree  of 
refraction  is  the  refractive  index  (n).  A  time  averaged,  or  "standard," 
atmosphere  will  bend  a  horizontally  propagating  electromagnetic  wave  downward 
with  a  curvature  less  than  the  curvature  of  the  earth's  surface.  If  the  re¬ 
fraction  is  sufficiently  great  to  cause  a  downward  curvature  that  exceeds  the 
earth's  surface  curvature,  then  the  propagating  electromagnetic  wave  will  be 
channeled  in  a  "duct."  Close  to  the  earth's  water  surface,  less  than  approx¬ 
imately  30  meters,  atmospheric  refraction  is  dominated  by  evaporation  and  any 
duct  so  formed  is  called  an  evaporation  duct. 

Typically,  atmospheric  water  vapor  content  over  ocean  areas  will  rapidly 
decrease  with  height.  As  a  result,  at  some  height,  the  refractive  index  gra¬ 
dient  will  decrease  below  that  value  required  for  ducting.  The  height  of  this 
critical  refractive  index  gradient  value  is  known  as  the  evaporation  duct 
height  (6). 

Evaporation  ducts  are  routinely  observed  over  the  world's  oceans,  with 
increasing  frequency  and  larger  duct  heights  being  found  at  low  latitudes, 
during  the  summer  season,  and  during  daylight  hours. 

Since  the  channeling  of  electromagnetic  energy  can  lead  to  greatly  ex¬ 
tended  ranges  for  radar  detection  of  surface  targets,  and  conversely,  greater 
counter-detection  ranges,  knowledge  of  n  and  5  •  becomes  vital  for  tactical 
operations.  In  addition,  6  serves  as  an  input  parameter  to  various  mathemati¬ 
cal  models  to  predict  pathloss,  and  thus,  a  detection/communication  range. 

In  practice,  direct  measurement  of  6  is  not  possible  with  a  radiosonde  or 
refractometer.  The  duct  height  must  be  inferred  from  meteorological  processes 
occurring  at  the  air/ocean  interface. 

PURPOSE 

The  purpose  of  the  work  reported  here  was  to  accomplish  four  objectives. 
Objective  I 

Evaluate  the  current  meteorological  measurement  technique  to  determine  if 
it  is  sufficient  to  adequately  describe  the  surface  meteorological  processes 
used  to  infer  6 . 

Objective  II 

To  compare  relative  performance,  sensitivities  to  meteorological  inputs, 
and  ease  of  computation  for  the  evaporation  duct  height  models  developed  by 
Jeske  (ref  1),  Ro the ran  (ref  2),  and  Davidson  (ref  3). 


Objective  ill 

To  compare  and  evaluate  the  EM  wave  propagation  pathless  models  developed 
by  Rotheram  (ref  2)  and  the  Naval  Ocean  Systems  Center  (NOSC)  (ref  4). 

Objective  IV 

To  determine  a  maximum  range  of  detection  error  for  the  modified  NOSC 
propagation  model  as  employed  by  the  Integrated  Refractive  Effects  Prediction 
System  (IREPS). 


OBJECTIVE  I 


The  theoretical  framework  for  calculation  of  8  used  by  all  three  models 
under  consideration  is  proposed  by  Monin  and  Obukhov  (ref  5)  and  is  based  upon 
the  relation  of  surface  layer  profiles  of  temperature,  moisture,  and  turbu¬ 
lence  with  surface  fluxes  of  mosientum,  sensible  heat,  and  latent  heat. 

Multilevel  measurements  needed  to  determine  these  surface  fluxes  are  ex¬ 
tremely  difficult  and  are  operationally  beyond  the  capability  of  field  meteor¬ 
ologists  around  the  world.  For  this  reason,  a  "bulk"  measurement  method  is 
sought.  If  the  Monin-Obukhov  expressions  for  temperature,  moisture,  etc.  are 
expressed  as  mathematical  derivatives,  then  the  integral  from  a  rough  surface 
boundary  (zq)  to  some  reference  height  (s)  will  express  the  fluxes.  ibis 

introduces  two  integration  boundary  conditions.  By  making  certain  assumptions 
about  the  meteorological  conditions  at  zQ,  only  the  four  variables  of  air 

temperature  (T  ),  sea-surface  temperature  (T  ),  moisture  (Rh) ,  and  wind  veloc- 
A  8 

ity  (U),  all  measured  at  z,  and  the  determination  of  z  need  to  be  made  to 

o 

determine  these  fluxes.  these  meteorological  measurements  then  are  referred 
to  as  "bulk"  measurements. 

The  models  developed  by  Davidson  uses  these  "bulk"  meteorological  measure¬ 
ments  to  explicitly  calculate  drag  coefficients  which  are  used  in  turn  to 
calculate  an  ata»spheric  stability  factor.  this  calculated  stability  factor 
is  used  to  define  a  stability  function.  From  the  stability  function,  the 
fluxes  of  temperature,  etc.  are  calculated.  Finally,  these  fluxes  are  used  in 
an  implicit  manner  to  calculate  8 . 

The  models  employed  by  Jeske  and  Ro the ram  uses  the  "bulk"  meteorological 

measurements  to  calculate  a  refractive  difference  between  z  and  z  ,  a  scale 

o 

length  empirically  derived  from  the  Bulk  Richard's  number,  and  a  profile  coef¬ 
ficient.  Based  upon  two  empirically  derived  stability  functions,  the  refrac¬ 
tive  difference  and  the  scale  length  are  used  to  explicitly  calculate  8 . 

To  achieve  the  first  objective  of  evaluating  the  "bulk"  meteorological 
measurements  technique,  observational  meteorological  and  pathloss  data  gath¬ 
ered  during  an  extensive  experimental  program  in  the  eastern  Mediterranean 
(ref  6)  are  used.  During  the  four  months  of  February,  April,  August,  and 
Movember  of  1972,  radio  propagation  measurements  for  frequencies  of  1.0426  GHz 
(L-band),  3.0075  GHz  (S-band),  9.624  GHz  (X-band),  17.9648  GHz  (KU-band),  and 
37.44  GHz  (KA-band)  were  taken  across  the  Aegean  Sea  between  the  islands  of 
Naxos  and  Mykonos,  a  horizontal  path  of  35.2  km.  For  a  transmitter  height  of 
approximately  5  meters  above  mean  sea  level  (MSL)  and  receiver  heights  of 
approximately  5,  10,  and  19  meters  above  MSL,  propagation  measurements  were 
made  at  15-minute  intervals. 

Meteorological  measurements  of  T#,  T#,  Rh,  and  U  were  collected  a  hourly, 
2-hourly  and  6-hourly  intervals,  depending  on  the  month,  for  transmitter 


and  receiver  sites.  In  addition,  meteorological  measurements  along  the  prop¬ 
agation  path  were  recorded  from  a  small  fishing  boat. 

As  discussed  earlier,  the  complex  processes  occurring  at  the  air/ ocean 
boundary  demand  sophisticated  measurements  for  detailed  description.  Since 
the  neeeesary  degree  of  sophistication  is  beyond  the  capability  of  field  mete¬ 
orologists  simple  "bulk"  meteorological  observations  were  purposely  made. 

Hie  evaporation  duct  height  models  of  Jeske,  Ro the ram,  and  Davidson  are 
initiated  with  these  data.  The  resultant  5  for  each  model  is  then  associated 
with  a  propagation  pathloss  measured  directly  and  at  the  same  time  as  the 
meteorological  measurement.  The  duct  height  and  its  associated  pathloss  are 
then  compared  to  a  theoretical  distribution  as  described  by  reference  4. 

OBJECTIVE  II 

For  the  second  objective,  the  "bulk" .meteorological  observations  are  used 
to  drive  each  model  to  produce  a  plot  of  5  with  time.  Differences  in  computed 
6  between  the  three  models  are  then  compared  for  relative  performance  and 
sensitivity  to  "bulk"  meteorological  parameters. 

OBJECTIVE  III 

To  achieve  the  third  objective,  the  propagation  pathloss  models  of  Rothe- 
ram  and  NOSC,  hereafter  referred  to  as  the  IREPS  model,  are  employed. 

Propagation  pathloss  calculations  from  the  IREPS  model  are  conducted  in 
two  steps.  First,  the  Jeske  evaporation  duct  height  model  is  employed  for 
determination  of  an  n  profile  under  neutral  conditions.  Neutral  conditions 
are  defined  as  those  where  the  air-sea  surface  temperature  difference  is  equal 
to  zero.  Secondly,  using  this  piecewise  linear  n  profile,  numerical  tech¬ 
niques  are  employed  to  calculate  reflection  coefficients.  Numerical  solutions 
of  the  fundamental  equation  of  mode  theory  are  then  found  using  these  reflec¬ 
tion  coefficients. 

The  Rotheram  model  for  propagation  loss  also  seeks  to  solve  the  funda¬ 
mental  equation  of  mode  theory. 

Differences  of  concern  between  the  two  models  come  from  two  factors. 

1 .  Hie  IREPS  model  assumes  no  effect  from  atmospheric  stability  on  the 
propagation  loss.  Rotheram1 s  model,  on  the  other  hand,  includes  a  stability 
dependent  effective  earth's  radius  factor  in  the  calculation  of  the  standard 
radio  horizon. 

2.  Hie  IREPS  model  assumes  a  smooth  earth  surface,  whereas  the  Rotheram 
model  includes  a  root-mean-square  sea  roughness  height  (standard  deviation  of 
the  sea  surface  from  its  mean  level)  in  the  computation  of  the  attenuation 
value . 

Each  model  is  driven  with  the  "bulk"  meteorological  observation  as  dis¬ 
cussed  earlier.  The  resultant  propagation  loss  is  compared  to  the  observed 
propagation  loss  associated  with  the  meteorological  observation.  A  statis¬ 
tical  distribution  of  propagation  pathloss  also  is  created  for  each  model  and 


com pared  with  the  statistical  distribution  of  observed  pathloss  values.  This 
■ smoothing"  allows  for  comparison  of  theory  without  undue  influence  of  "ex¬ 
treme*  input  parameters. 

OBJECTIVE  XV 

To  achieve  the  fourth  objective,  the  IREPS  model  is  driven  with  the 
"bulk”  meteorological  observations.  An  estimated  pathloss  value  is  generated 
for  a  fixed  range  of  35.2  km,  the  range  between  the  transmitter  and  receiver 
for  pathloss  observations.  This  estimated  pathloss  is  in  turn  used  to  calcu¬ 
late  an  estimated  maximum  detection  range.  By  using  the  observed  pathloss 
value,  an  "inferred"  actual  maximum  detection  range  is  calculated.  A  compar¬ 
ison  of  estimated  maximum  range  and  inferred  maximum  range  is  made  to  deter¬ 
mine  detection  range  error  boundaries  for  use  in  assessing  IREPS  reliability. 
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RESULTS 


CURRENT  TECHNIQUE  EVALUATION 

The  first  two  objectives  of  this  report  are  met  by  considering  the  evapo¬ 
ration  duct  height  models  of  Jeske,  Ro the ram,  and  Davidson  and  employing  the 
meteorological  observations  for  all  four  time  periods  and  both  island  loca¬ 
tions.  Because  of  the  frequency  of  November  meteorological  data  and  because 
relative  model  performance  is  similar  for  each  period,  the  November  data  will 
be  used  in  all  comparisons  unless  stated  otherwise.  Also,  for  comparative 
purposes,  a  reasonable  5  will  be  defined  as  40  meters  or  less. 

Figures  1  and  3*  illustrate  6  as  a  function  of  time  for  each  model.  Gen¬ 
eral  trends  of  synoptic  and  diurnal  height  fluctuations  are  nearly  equally 
produced  by  all  three  models,  but  performance  on  individual  data  shows  a  great 
variance  and  will  account,  in  some  part,  for  the  extremes  in  comparison  data. 
For  example,  identical  meteorological  inputs  produced  5  of  16,  30,  and  1189 
meters  by  the  models  of  Jeske,  Rotheram,  and  Davidson,  respectively,  ttie  dis¬ 
continuity  in  the  Davidson  5  data  is  due  to  the  model's  inability  to  handle 
certain  meteorological  combinations  and  will  be  explained  later.  Table  1 
displays  the  basic  statistics  for  the  evaporation  duct  height  as  calculated  by 
the  models  of  Jeske,  Rotheram  and  Davidson. 

Davidson's  Model 

Of  the  three  evaporation  duct  height  models,  Davidson's  is  the  most  vari¬ 
able.  Duct  heights  ranged  between  0  and  1189  meters  with  a  mean  6  of  22.7 
meters.  Correlation  with  the  5  of  Jeske  and  Rotheram  is  0.273  and  0.365, 
respectively. 

For  stable  conditions,  air  temperature  warmer  than  sea  surface  tempera¬ 
ture,  an  air-sea  temperature  difference  of  1 °C  coupled  with  low  humidities  and 
wind  velocities,  is  sufficient  to  produce  an  unrealistic  6  .  For  example,  for 

Rh  *  55%  and  U  =  1  m  sec  1 ,  6  =  242  meters.  Increasing  Rh  by  10%  or  U  by  1  m 

sec  ,  decreases  6  to  14  and  49  meters,  respectively. 

For  neutral  conditions  and  U  *  1  m  sec  1 ,  the  model  failed  to  produce  a 
5  as  a  function  of  Rh  and  T  .  For  example,  if  T  =  T  =  19°C  and  Rh  *=  62%, 

3  21  S 

the  model  fails.  If  T  =  T  =  15°C,  the  model  will  fail  when  Rh  =  51%.  For  U 
-1  as 

>  1  m  sec  ,  the  model  produces  realistic  6  regardless  of  Rh. 

For  unstable  conditions,  air  temperature  colder  than  sea  surface  tempera- 
tures,  and  U  =  1  m  sec  ,  the  model  fails  to  produce  a  6  regardless  of  Rh.  If 
U  >  1  m  sec  ,  however,  reasonable  duct  heights  were  obtained  in  all  cases. 

For  the  November  observation  period,  the  meteorological  conditions  neces¬ 
sary  to  produce  an  unrealistic  6  or  a  model  failure  occurred  3.3  and  14.0 
percent  of  the  time,  respectively.  For  the  February  period,  the  percentages 


•All  figures  and  tables  are  placed  at  the  end  of  the  report. 
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ara  5.1  and  13.6;  for  the  April  period,  the  percentages  are  21.1  and  25.0;  and 
for  the  August  period,  unrealistic  5  or  failure  percentages  are  26.7  and  14.0, 
respectively. 

Rotheram'*  Model 

Rotheram' s  evaporation  duct  height  model  is  much  less  variable  with  in¬ 
dividual  meteorological  data  than  that  of  Davidson.  £  varies  from  0  to  253 
meters  with  a  mean  £  of  11.5  meters.  £  correlations  with  Jeske's  and  David¬ 
son's  models  are  0.93  and  0.365,  respectively. 

As  with  Davidson's  model,  stable  conditions  coupled  with  low  humidities 
and  wind  velocities  are  sufficient  to  produce  an  unrealistic  £  ,  but  not  to  the 
extent  shown  by  the  Davidson  model.  For  example,  with  T  =  20°C,  T  =  19°C, 

^  3  8 

Rh  -  78%,  and  U  *  3  m  sec  ,  6-182  and  327  meters  for  Rotheram  and  Davidson, 
respectively. 

For  neutral  and  unstable  conditions,  the  model  never  failed  to  produce  a 
£,  regardless  of  the  meteorological  parameters . 

The  meteorological  conditions  necessary  to  produce  an  unrealistic  fi  oc¬ 
curred  8.4,  5.7,  4.9,  and  3.9  percent  of  the  time  for  February,  April,  August, 
and  November,  respectively. 

Jeske's  Model 

Jeske's  evaporation  duct  height  model  was  the  least  sensitive  to  individ¬ 
ual  data,  fi  varies  between  0  and  134  meters  with  a  mean  6  of  10.7  meters. 

In  all  stability  cases,  Jeske's  model  fails  to  produce  a  fi  only  when  u  = 
0.  This  occurred  3.4,  5.7,  0,  and  0  percent  of  the  time  for  February,  April, 
August,  and  November,  respectively. 

The  Jeske  model  behaves  very  similar  to  that  of  Rotheram.  For  stable 

conditions,  low  humidities  and  wind  velocities  produce  unrealistic  fi  ,  but  not 

to  the  extent  shown  by  Rotheram.  For  example,  if  T  -  20°C,  T  »  19°C,  Rh  « 

as 

78%,  and  U  -  3  n  sec  ,  the  6-80  and  1 82  meters  for  Jeske  and  Rotheram, 
respectively. 

Meteorological  conditions  necessary  to  produce  an  unrealistic  5  occurred 
11.8,  38.4,  52.1,  and  2.0  percent  of  the  time  for  February,  April,  August,  and 
November,  respectively. 

For  neutral  and  unstable  conditions,  reasonable  fi  were  obtained  100  per¬ 
cent  of  the  tine  regardless  of  the  meteorological  parameters,  except  for  0  * 
0,  as  previously  discussed. 

PERFORMANCE,  METEOROLOGICAL  INPUT,  AMD  COMPUTATION  COMPARISONS 

Because  of  the  extreme  complexity  of  measuring  6  directly,  actual  5  val¬ 
ues  ware  not  obtained  during  the  observation  period.  Therefore,  the  duct 
heights  calculated  by  each  model  are  compared  to  a  theoretical  distribution  of 


8  as  a  function  of  pathloss,  as  discussed  in  reference  4.  Figures  4  through 
IS  show  the  distribution  of  calculated  6  versus  observed  pathloss  at  four 
frequencies  for  each  duct  height  model.  Superimposed  on  each  diagram  is  the 
theoretical  dependence  of  5  with  pathloss. 

From  visual  inspection,  there  appears  to  be  no  significant  difference  in 
the  performance  of  the  three  models.  While  there  is  general  agreement  with 
the  theoretical  curves  for  each  frequency,  there  exists  a  great  variability 
within  the  individual  data  points,  with  the  greatest  variability  occurring  at 
low  duct  heights. 

A  mean  pathloss  error  is  created  by  taking  the  mean  difference  between 
observed  and  theoretical  pathloss  for  each  calculated  duct  height.  The  mean 
pathloss  error  as  a  function  of  6  and  frequency  is  shown  in  figures  16  through 
19.  Again,  no  significant  difference  exists  between  the  duct  height  models. 
Compared  to  theoretical  pathloss  values,  there  appears  to  be  less  than  or 
equal  to  a  10-dB  error  for  each  model  except  for  ducts  less  than  7  meters, 
where  all  three  models  show  an  inferred  overestimated  pathloss. 

While  all  three  models  show  a  mean  statistical  agreement  with  present 
theory,  model  sensitivities  to  combinations  of  meteorological  parameters  pro¬ 
duce  an  extreme  variation  of  5  .  The  assumption  of  "bulk"  meteorological  param¬ 
eters  being  sufficient  for  accurate  individual  6  determination  is  called  into 
suspicion.  It  is  seen  that  a  Tg  variation  of  just  1  °C  can  produce  totally 

different  duct  heights.  Under  normal  operations,  meteorological  personnel  are 
unable  to  measure  the  sea  surface  temperature  to  an  accuracy  of  l*C.  In  addi¬ 
tion,  given  the  influences  of  a  ship's  superstructure  and  operations  such  as 
launching  aircraft,  the  other  parameters  of  T  ,  Rh,  and  u  are  equally 

difficult  to  measure  to  an  accuracy  required  by  any  model. 

The  needed  accuracy  of  meteorological  inputs  also  calls  into  question  the 
assumption  of  horizontal  atmospheric  homogeneity.  Any  synoptic  or  mesoscale 
circulation  that  would  lead  to  horizontal  homogeneity  would  also  tend  to  pro¬ 
duce  neutral  or  stable  conditions.  It  is  seen  that  all  models  display  their 
greatest  sensitivities  within  the  stable  regimes.  For  unstable  regimes,  it  is 
not  unrealistic  to  observe  both  temporal  and  spatial  meteorological  parameter 
fluctuations  greater  than  that  required  by  the  models.  The  pathloss  errors 
experienced  between  observed  and  inferred  pathlosses  clearly  reflect  a  desir¬ 
ability  to  consider  atmospheric  inhomogeneities  over  the  propagation  path. 
Current  operational  propagation  models  are  incapable  of  addressing  this  topic, 
however.  An  extremely  expensive  and  difficult  research  and  development  effort 
would  be  needed  to  collect  and  analyze  data  and  to  develop  new  propagation 
algorithms.  Also,  it  is  operationally  impractical  to  sample  the  atmosphere  in 
i  detail  sufficient  to  adequately  describe  the  propagation  path. 

COMPARISON  AND  EVALUATION  OF  PATHLOSS  MODELS 

The  third  objective  of  this  report  is  to  compare  and  evaluate  the  propa¬ 
gation  loss  models  of  Rotheram  and  IREPS. 

Figures  20  through  29  are  scatter  plots  of  calculated  pathloss  from  the 
models  of  Rotheram  and  IREPS  versus  the  observed  pathloss  for  the  November 
observation  period.  Tables  2  and  3  display  the  basic  statistics  for  path- 
losses  as  observed  and  calculated  by  the  pathloss  models  of  IREPS  and  Rotheram. 
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At  the  L-band  frequency,  the  standard  deviations  of  observed,  I REPS,  and 
Rotheraa  pathloss  are  3.5  dB,  4*0  dB,  and  11.0  dB,  respectively.  As  the  fre¬ 
quency  increases,  the  standard  deviation  of  all  three  path losses  increases, 
with  the  Rotheraa  standard  deviation  increasing  most  rapidly.  At  the  KA-band 
frequency,  the  observed,  I REPS,  and  Rotheraa  pathloss  standard  deviations  are 
5.6  dB,  15.7  dB,  and  24.8  dB,  respectively. 

The  correlation  between  the  IREPS  and  Rotheraa  pathloss  aodels  is  high 
for  frequencies  of  x-band  and  below.  For  exaaple,  at  S-band,  the  correlation 
between  pathlosses  is  0.93.  At  frequencies  above  X-band,  however,  the  correla¬ 
tion  falls  to  zero,  ttiis  can  be  explained  by  the  ZREPS  aodel  not  considering 
surface  roughness,  which  beccaes  increasingly  inportant  with  high  frequencies. 

The  correlation  between  both  aodel  pathlosses  and  the  observed  pathloss 
increases  froa  0.35  at  L-band  to  0.59  at  X-band,  but  again  falls  to  zero  at 
KA-band  frequency.  This  is  expected,  however,  since  both  models  were  devel¬ 
oped  to  operate  below  20  GHz,  thereby  eliminating  consideration  of  atmospheric 
EM  wave  absorption. 

Figures  30  through  34  show  the  cumulative  percent  occurrence  of  observed 
pathloss  with  the  pathlosses  calculated  by  Rotheram  and  IREPS  as  a  function  of 
frequency.  The  general  trends  of  the  curves  taken  together  with  the  statisti¬ 
cal  summary  demonstrate  the  statistical  agreement  with  theory,  but  as  with  the 
evaporation  duct  height  models,  a  large  variance  in  individual  data  leads  to 
discrepancies  between  model  performance  and  observational  data.  Investigation 
of  time  periods  other  than  November  show  similar  statistical  trends  and  model 
variabilities. 

Statistically,  the  IREPS  modeled  pathloss  is  less  variable  than  the  Rothe¬ 
raa  modeled  pathloss,  and  therefore  more  closely  represents  the  pathloss  dis¬ 
tribution  as  observed.  Since  a  major  input  to  the  model  is  6  ,  and  Jeske's 
duct  height  model  employed  by  the  IREPS  pathloss  model  was  the  least  variable, 
this  is  expected. 

The  correlation  of  the  IREPS  calculated  pathloss  to  observed  pathloss 
reaches  a  maximum  at  approximately  10  GHz,  but  is  poor  above  this  frequency. 
At  higher  frequencies,  where  surface  roughness  should  influence  the  pathloss 
models,  the  correlation  between  both  IREPS  and  Rotheram' s  pathloss  models 
output  and  observed  pathloss  is  not  significantly  sufficient  to  justify  inclu¬ 
sion  of  a  surface  roughness  parameter  as  defined  by  Rotheraa. 

MAXIMUM  DETECTION  ERROR  DETERMINATION 

The  final  abjective  of  this  report  is  to  determine  a  maximum  radar  detec¬ 
tion  range  prediction  error  for  the  IREPS  propagation  model. 

Figures  35  through  39  are  scatter  plots  for  the  estimated  maximum  detec¬ 
tion  range  calculated  by  the  IREPS  model  versus  an  inferred  maximum  detection 
range  based  upon  an  observed  pathloss  over  the  measurement  range.  Both  esti¬ 
mated  and  inferred  maximum  detection  ranges  are  based  on  a  maximum  free  space 
detection  range  of  185  km. 

Both  calculated  and  inferred  maximum  detection  ranges  vary  greatly  be¬ 
cause  of  the  previously  demonstrated  sensitivity  of  5  to  meteorological  param- 


eters,  and  the  sensitivity  of  the  propagation  models  height-gain  function 
relationship  to  5  .  Figure  40  is  an  example  of  the  height-gain  maximization 
curves  used  to  extrapolate  an  inferred  detection  range  from  an  observed  path- 
loss,  and  illustrates  the  height-gain  function,  6  relationship.  Xt  becomes 
evident  from  the  figure  that  a  very  small  pathloss  change,  due  to  a  change  in 
6,  can  cause  an  extreme  change  in  maximum  inferred  detection  range. 

Correlations  between  calculated  and  inferred  maximum  detection  ranges 
varied  from  a  -0.02  at  3  GHz  to  a  maximum  of  0.5  at  10  GHz.  By  calculating 
the  difference  between  computed  and  inferred  maximum  detection  ranges,  divid¬ 
ing  this  difference  by  the  computed  range,  and  then  segregating  these  relative 
range  errors  into  percent-error  categories,  a  histogram  of  percent  occurrence 
versus  relative  detection  range  error  is  derived.  The  histograms  for  each 
frequency  band  are  displayed  in  figures  41  through  45.  For  L-  and  S-band 
frequencies  that  are  not  significantly  affected  by  the  evaporation  duct,  rel¬ 
ative  detection  range  errors  of  less  than  50%  occurred  98.1%  and  90.7%  of  the 
time,  respectively.  For  X-band  frequencies,  the  relative  detection  range 
errors  of  less  than  50%  occur  57.0%  of  the  time.  At  K-band,  where  the  evapo¬ 
ration  duct  has  the  greatest  influence  in  propagation,  the  IREPS  model  pro¬ 
duces  detection  range  errors  of  less  than  50%  only  30.8%  of  the  time.  Rel¬ 
ative  detection  range  errors  were  in  excess  of  100%  45.7  percent  of  the  time. 
For  the  KA-band  frequencies,  the  relative  detection  range  errors  of  less  than 
50%  increases  to  62.4%  of  the  time.  This  increase  is  to  be  expected,  since 
once  again  the  evaporation  duct  effect  diminishes.  The  increase  is  not  as 
dramatic  as  at  lower  frequencies,  however,  since  atmospheric  absorption,  which 
is  not  modeled  in  the  IREPS  program,  becomes  a  substantial  factor  in  propa¬ 
gation.  with  the  exception  of  the  KA-band  frequency,  the  IREPS  model  tends  to 
underestimate  a  maximum  detection  range,  with  the  greatest  underestimation 
occurring  at  S-band  frequencies. 

Statistical  analysis  of  the  XREPS  detection  range  software  cannot  be  used 
to  support  the  soundness  of  propagation  theory  as  well  as  it  did  evaporation 
duct  theory.  The  model's  sensitivity  to  the  evaporation  duct  height  is  viv¬ 
idly  demonstrated  by  the  frequency  dependent  detection  range  errors,  but  the 
ability  of  the  model  to  perform  when  the  influence  of  the  evaporation  duct 
height  is  greatest  is  poor. 

For  this  reason,  it  becomes  even  more  important  to  understand  that  the 
IREPS  software  is  written  to  assess  relative  emitter  systems  performance,  and 
not  to  produce  an  absolute  propagation  range  value. 


CONCLUSIONS 


t.  The  evaporation  duct  height  aodels  of  Jeske,  Rotheram,  and  Davidson  func¬ 
tion  equally  well  with  reasonable  duct  heights  calculated  when  considering  the 
statistical  average  of  meteorological  inputs.  When  individual  meteorological 
data  points  are  considered,  however,  the  three  models  show  a  wide  fluctuation 
in  outputs.  The  model  of  Jeske  shows  the  most  reasonable  duct  height  for  any 
particular  meteorological  situation,  while  the  model  of  Davidson  shows  the 
greatest  variance.  Because  of  the  equal  functioning  when  considering  statis¬ 
tical  average,  the  model  performance  is  more  a  function  of  the  individual 
meteorological  parameters  than  of  any  one  technique  of  duct  height  calculation. 

2.  The  meteorological  observations  taken  by  operational  personnel  are  not  of 
sufficient  quality  to  insure  proper  duct  height  calculations  when  individual 
duct  height  predictions  are  desired. 

3.  The  evaporation  duct  height  is  a  major  input  to  the  propagation  models  of 
IREPS  and  Rotheram.  Since  individual  duct  heights  show  large  variances,  both 
models  reflect  large  fluctuations  in  individual  pathloss  calculations.  When 
considering  statistical  means  in  observational  data,  both  models  function 
reasonably  well,  with  the  IREPS  model  more  closely  representing  observation. 
Since  the  IREPS  model  does  not  include  a  direct  consideration  of  surface  rough¬ 
ness  as  does  Rotheram' s  model,  but  performs  equally  or  slightly  better  within 
same  frequency  bands,  the  additional  effort  to  model  surface  roughness  does 
not  appear  to  be  justified. 

4.  A  second  consequence  of  the  model's  sensitivity  to  individual  meteoro¬ 
logical  inputs  involves  a  reasonable  model  performance  when  assuming  atmo¬ 
spheric  horizontal  hcaiogeneity.  It  is  obvious  that  in  regions  of  air  mass 
boundaries,  horizontal  homogeneity  is  a  poor  assumption.  It  is  shown  that  the 
models'  sensitivities  to  meteorological  inputs  indicates  a  desirability  for 
consideration  of  a  heterogeneous  atmosphere .  Given  the  inability  of  meteoro¬ 
logical  personnel  to  make  sufficiently  accurate  measurements,  however,  at¬ 
tempts  to  develop  new  theoretical  techniques  to  allow  more  inputs  for  better 
horizontal  resolution  would  prove  counterproductive.  Assuming  model  sensi¬ 
tivities  to  meteorological  input  could  be  overcame  and  new  propagation  algo¬ 
rithms  for  a  heterogeneous  atmosphere  could  be  developed,  the  practicality  and 
expense  of  operationally  gathering  sufficient  quantities  of  input  data  to 
significantly  reflect  atmospheric  conditions  would  be  unreasonable.  When 
considering  all  problems  involved,  the  assumption  of  a  homogeneous  atmosphere 
for  a  fleet  operational  propagation  model  is  realistic. 

5.  All  previous  conclusions  indicate  the  theory  behind  calculation  of  duct 
height  and  pathloss  values  is  statistically  sound.  For  statistical  studies  of 
propagation  phenomena  or  development  of  electromagnetic  systems,  both  the 
XREPS  and  Rotheram  models  perform  reasonably  well  and  are  accurate  predictors. 
However,  the  use  of  the  XREPS  model  for  estimation  of  individual  propagation 
ranges  based  upon  an  individual  meteorological  observation  is  suspect.  At  L-, 
S-,  X-,  KU-,  and  KA-band  frequencies,  the  IREPS  model  produced  median  relative 
range  errors  of  5,  19,  40,  83,  and  31  percent,  respectively.  The  poorest  per¬ 
formance  of  IREPS  occurs  at  frequencies  most  affected  by  the  evaporation  duct. 
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RECOMMENDATIONS 


1 •  Qiality  of  individual  meteorological  observations  needed  as  input  to 
current  state-of-the-art  propagation  models  is  insufficient.  the  quality  can 
be  improved  by  two  methods: 

a.  Improve  the  observational  skills  of  the  individual  through  further 
training,  with  emphasis  on  attention  to  established  observational  procedures. 

b.  Develop  a  sensor  system  that  can  be  deployed  away  from  heating/cool¬ 
ing  and  turbulent  effects  created  by  the  observation  platform. 

2.  In  the  face  of  insufficient  quality  in  observed  meteorological  data  and 
the  sensitivity  of  duct  height  models  to  these  data,  develop  an  entirely  new 
duct  height  model  using  statistical  techniques. 

3.  Current  propagation  models  are  based  upon  pathloss  values  created  from 
assumed  conditions  of  atmospheric  absorption,  atmosheric  homogeneity,  system 
noise,  etc.  Controlled  experiments  need  to  be  designed  to  investigate  the 
role  of  each  assumption  upon  the  actual  pathloss  over  the  propagation  path. 

4.  Train  fleet  operators  of  the  I REPS  model  in  the  limitations  and  recog¬ 
nition  of  situations  in  which  the  model's  performance  is  marginal  or  poor. 

5.  Undertake  further  evaluation  of  the  IREPS  model  by  comparing  a  predicted 
range  with  an  actual  observed  maximum  detection  range. 
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Table  1 .  Basic  statistics  for  the  evaporation  duct  height  as  calculated  by 
the  models  of  Jeske,  Rotheram,  and  Davidson 


Source 


Standard  Deviation 
(m) 


Maximum 

(m) 


Jeske  model 

10.7 

13.4 

134.5 

Rotheram  model 

11.5 

23.7 

353.7 

David  so  n  model 

22.7 

96.8 

1189.9 

Correlation  Matrix 

Jeske  model  Rotheram  model 

Davidson  Model  0.2729580 

0.3655519 

Jeske  model 

0.9281335 

Minimum 

(m) 


Table  3.  Correlation  of  pathlosses  as  observed  and  calculated  by  the  pathloss 
model 8  of  XREPS  and  Rotheram. 


Frequency 


Correlation  Matrix 


Observed 


Rotheram 

model 


.624  GHz  Observed  Path  Loss  (dB) 
Receive  Antenna  Height  *  10  (m) 


17.9648  GHz  Observed  Path  Loss  (dB) 
Receive  Antenna  Height  »  9.5  (m) 


Receive  flntenn 


Duct  Height  (m).  Calculations  Based  On 
NPS6 1-78-005  of  15  Rug  78  (Davidson) 

Figure  1 2.  Davidson  evaporation  duct  height  versus  L-band  frequency  observed  pathloss. 


Mykonos  Meteorological  Data  -  November  1972 


Pathloss  theory  curve  from  NELC/TR  194 


Duct  Height  (m).  Calculations  Based  On 
NPS6 1-78-005  of  15  Aug  78  (Davidson) 

Figure  13.  Davidson  evaporation  duct  height  versus  S-band  frequency  observed  pathloss. 


17.9648  GHz  Observed  Path  Loss  (dB) 
Receive  Antenna  Height  *  9.5  (m) 


Mykonos  Meteorological  Data  -  November  1972 


Pathloss  theory  curve  from  NELC/TR  194 


Duct  Height  (m).  Calculations  Based  On 
NPS6 1 -78-005  of  15  Aug  78  (Davidson) 

Figure  15.  Davidson  evaporation  duct  height  versus  KU-band  frequency  observed  pathloss. 


Mykonos  Meteorological  Data  -  November  1972 


Model  underestimated  pathloss 


Model  overestimated  pathloss 


Duct  Height  (m).  Calculations  Based  On  Theory  Of 
Jeske  (solid),  Rotheram  (dot),  and  Davidson  (dash) 

Figure  16.  Three  model  evaporation  duct  height  versus  L- band  frequency  mean  pathloss  error. 


Mykonos  Meteorological  Data  -  November  1972 


Model  underestimated  pathloss 


Model  overestimated  pathloss 


Duct  Height  (m).  Calculations  Based  On  Theory  Of 
Jeske  (solid),  Rotheram  (dot),  and  Davidson  (dash) 

Figure  18.  Three  model  evaporation  duct  height  versus  X-band  frequency  mean  pathloss  error. 


GHz  Obs< 
ve  Rnter 


GHz  Obse 
i  ve  flnter 


17.9648  GHz  Calculated  Path  Loss  CdB. ) 

Based  on  NOSC/TR  771  of  28  Rpr  82 

Figure  23.  IREPS  KU-band  frequency  pathloss  versus  observed  KU-band  frequency  pathloss. 


Mykonos  Meteorological  Data  -  November  1972 
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.0426  GHz  Obs< 
Receive  flnter 


GHz  Observed  Path  Loss  (dB) 


9.624  GHz  Calculated  Path  Loss  (dB) 

Based  on  Marconi  Lab  MIR. 77/37 

Figure  27.  Rotheram  X-band  frequency  pathloss  versus  observed  X-band  frequency  pathloss 


GHz  Obse 
ve  flnten 


Figure  31 .  S-band  frequency  pathloss  versus  cumulative  percent  occurrence. 
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Predicted  Range  (Km) 

Based  on  NOSC/TR  771  of  28  Rpr  82 

Figure  35.  IREPS  predicted  range  versus  inferred  actual  range  for  L-band  frequency 


Inferred  fl( 
Based  on  Ob 


3.0075  GHz  -  Receiver  Height  10  M 


Figure  36.  IREPS  predicted  range  versus  inferred  actual  range  for  S-band  frequency. 
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Predicted  Range  (Km) 

Based  on  NOSC/TR  771  of  28  Apr  82 

Figure  38.  IREPS  predicted  range  versus  inferred  actual  range  for  KU-band  frequency 


Predicted  Range  (Km) 

Based  on  NOSC/TR  771  of  28  Apr  82 

Figure  39.  IREPS  predicted  range  versus  inferred  actual  range  for  KA-band  frequency 
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Figure  40.  Observed  X-band  frequency  pathloss  versus  inferred  propagation 
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